Compression of Periodic Complementary 
Sequences and Applications!^ 

Dragomir Z. Dokovicj?], Ilias S. Kotsirea^ 

Abstract 

A collection of complex sequences of length v is complementary if the sum of their pe- 
riodic autocorrelation function values at all non-zero shifts is constant. For a complex 
^ . sequence A = [aQ,ai, . . . ,ay-i] of length v = dm we define the m-compressed sequence 

Qj [ A^^^ of length d whose terms are the sums + aj+rf + • • • + aj+(m_i)rf. We prove that 

PL^ ■ the m-compression of a complementary collection of sequences is also complementary. 

The compression procedure can be used to simplify the construction of complementary 
{±l}-sequences of composite length. In particular, we construct several supplementary 
Q I difference sets (SDS) {v;r,s;X) with v even and A = (r + s) — v/2, given here for the 

' first time. There are 15 normalized parameter sets {v;r,s;X) with v < 50 for which the 

existence question was open. We resolve all but one of these cases. 
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1 Introduction 

We consider the ring Z^, = {0, 1, . . . , f — 1} of integers modulo a positive integer v. Let ki, . . . ,kt 
be positive integers and A an integer such that 

t 

X{v-l) = J2Hki-l), (1) 

i=l 

and let Xi, . . . ,Xt be subsets of Zy such that 

\Xi\ = h, te{i,...,t}. (2) 

I Definition 1 We say that Xi, . . . , Xt are supplementary difference sets (SDS) with parameters 
{v; ki, . . . , kf, A), if for every nonzero element c G Z„ there are exactly A ordered pairs {a, h) such 
that a — b = c (mod v) and {a, b} C Xi for some i G {1,2, . . . ,t} . 

These SDS are defined over the cyclic group Zy. More generally SDS can be defined over 
any finite abelian group, and there are also further generalizations where the group may be any 
finite group. However, in this paper we shall consider only the cyclic case. 
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In the context of an SDS, say Xi, . . . , Xt, with parameters {v; ki, . . . , kt] A), we refer to the 
subsets Xi as the base blocks and we introduce an additional parameter, n, defined by: 

n= ki + --- + kt~ \. (3) 

The SDS formahsm is a generahzation of various well-known and studied designs. 

SDS consisting of just one base block, the case t = 1, are called cyclic difference sets and 
are denoted by (f ; k; A). The reader is referred to the classic book [2j, the more recent book 
[T7] and the recent survey [8]. The list of {v; k; A) cyclic difference sets with f < 50 in a normal 
form can be found in section 3 of [5j. 

The SDS with two base blocks, the case t = 2, covers several interesting designs. First, the 
circulant D-optimal designs [9l |6] corresponding to the case where v is odd and n = (t> — l)/2. 
Second, the binary complementary pairs (periodic analogs of complementary Golay pairs) [18] 
correspond to the case v = 2n. Third, when the two base blocks have the same size the SDS can 
be used to construct some BIBD (balanced incomplete block designs) [131 ITT] . We shall denote 
the SDS with two base blocks of size r and s by (f ; r, s; A). Upon imposing the normalization 
condition |>r>s>2, it can be seen that for f < 50 there are 227 feasible parameter sets, 
see [5]. These parameter sets are those that satisfy the conditions (jlj) and ([2]) with ki = r, 
k2 = s. There is a non-existence result for SDS {v;r,s;X) obtained in [1], which can be used 
to eliminate some of these 227 feasible parameter sets. In [12] the authors construct 30 new 
[v] r, s; A) and show that an additional 19 (f ; r, s; A) do not exist. In [5j the author constructs 
8 new (v;r, s; A). We list here the remaining 15 undecided cases from [5]. They all have v in 
the interval 40 < f < 50. 



(41;15,6;6) 


n= 


45 


(43;9,4;2) 


n= 


■11 


(44;19,2;8) 


n= 


■13 


(45;18,2;7) 


n= 


■13 


(46;21,6;10) 


n= 


■17 


(47;9,5;2) 


n= 


■12 


(47;12,3;3) 


n= 


■12 


(47;14,2;4) 


n= 


■12 


(47;15,5;5) 


n= 


■15 


(48;14,3;4) 


n= 


■13 


(49;10,3;2) 


n= 


■11 


(49;21,4;9) 


n= 


■16 


(50;8,7;2) 


n= 


■13 


(50;20,4;8) 


n= 


■16 


(50;22,21;18) 


n= 


■25 



In this paper (see Theorem [5]) we decide the existence of all of them except for the case 
(49; 21, 4; 9). The main tools that we use for that purpose are computational techniques based 
on the Power Spectral Density (PSD) test as well as the method of SDS compression. For 
the PSD test see section HI The compression method is studied systematically and in full 
generality (for arbitrary complex sequences and in particular for SDS) in section [5l We note 
that techniques that essentially amount to compression have been used previously [3|, [11] in the 
context of D-optimal designs. 

Apart from the cases where there exist Golay complementary pairs of length v, the binary 
complementary pairs are known to exist only for v = 34 [4j and = 50 [5], [lO] . In this paper we 
construct such pairs of length f = 50 and f = 58 (see the section [6]). The example for f = 50 
has different parameters, namely (50; 22, 21; 18), from the previously constructed examples with 
parameters (50; 25, 20; 20). The first unknown case is now for v = 68. 
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2 The group ring approach to SDS 



In the study of SDS it is often convenient to use the group ring of the additive group Z^, with 
coefficients in C, the field of complex numbers. This group ring can be identified with the 
quotient ring C[a;]/(x^ ~ 1)) where C[x] is the polynomial ring over C in a single indeterminate 
X and {x'" — 1) is the ideal generated by the polynomial x"" — 1. By abuse of notation, we shall 
consider x also as an element of this group ring, in which case we have x" = 1. Then Z^, is 
embedded in this group ring via the map which sends i — )■ a;*, i G Z^,. Note that because x^ = 1 
there is a ring homomorphism (evaluation at 1) C[x]/(x'' — 1) — > C which sends x — )■ 1. 

To any subset X C Z„ we shall associate an element of the group ring, namely X{x) : = 
Tlii&x abuse of notation, we shall denote this element also by X, except that in the case 

X = Zt, we set 

T = T{x) = 1 + X + ■■■ + x''-\ 

It will be clear from the context which meaning of X is used. 

The group ring has an involution which sends each complex number c to its complex con- 
jugate c and sends x to its inverse x~^ = x^~^. We shall denote this involution by an asterisk, 
e.g., we have x* = x~^ and T* = T. For any element X of the group ring we define its norm 
to be the element X(X) = XX*. We also define the function Nx : Z^, — )■ Z by declaring that 
Nx{s) is the coefficient of in N{X), i.e., we have 

v-1 

N{X) = J2Nx{t)x\ 

To any complex sequence of length v, say A = [oq, cti, . . . , we assign the element 
A{x) = Yli '^j^* of the group ring. Obviously, the map sending A to A{x) is injective. To 
simplify notation, we shall write A instead of A{x) if no confusion will arise. By evaluating 
A{x) at 1, we obtain A{1) = Y.i'^i- Similarly, we have N{A){1) = 

We point out that the SDS property is equivalent to an identity in the group ring. This is 
stated in the next lemma whose proof is straightforward. 

Lemma 1 Let Xi, . . . , Xt be subsets o/Z^ and assume that Eqs. ^\) and ^ hold. Then these 
subsets are the base blocks of an SDS with parameters {v; ki, . . . , kf, A) if and only if 

t 

J2N{X,)=n + XT, (4) 

i=l 

where n is defined by Eq. 

3 Power spectral density 

Many known SDS have been constructed by a computer search. In such searches it is often 
convenient to replace a subset X C Z„ with the {±l}-sequence A = [oq, ai, . . . , a^_i] where 
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tti = —1 if and only if i G X. We shall refer to A as the associated sequence of X. In that case 
we have 

A{x) =T{x) -2X{x). (5) 

Let A = [uq, ai, . . . , a„_i] be an arbitrary complex sequences of length v, and let u = e^'^*/" 
be a primitive v-th root of unity. 

The Discrete Fourier Transform (DFT) of the sequence A is the function Z„ — )■ C defined 
by the formula 

v-l 

DFT^(s) = ^a^w^". 

j=0 

We often identify a function on with the sequence of its values, e.g., we write 
DFT^ = B= [bo, 61, ... , 6.-1], where b, = DFTa(j)- 

This convention is used throughout the paper. 

The Power Spectral Density (PSD) of the same sequence A is the function — )■ R defined 
by the formula 

PSD^(s) = |DFT^(s)^ 
The Periodic Autocorrelation Function (PAF) of A is defined as 



v-l 
j=0 



We shall refer to the argument s as the shift variable. 

It is straightforward to verify that for any complex sequence A = [oq, ai, . . . , a„_i] and the 
corresponding element A = A{x) = Oq + aix + • • • + a^^ix^'^ of the group ring we have 

Na = PAFa. (6) 



In general we have Na{—s) = Na{s), and for a real sequence A we have Na{—s) = Na{s). 
Another important fact is the following classical theorem (see e.g. [1^ 



Theorem 1 (Wiener— Khinchin) For any complex sequence A = [qq, cti, . . . , Ot^-i], or equiv- 
alently any element A = + aix + ■ ■ ■ + a^^ix"""^ of the group ring, we have 

PSDa = DFT(PAF^). 

By using Eq. ([6]) this can be rewritten as follows: 

v-l 

PSD^(s) = ^Ar^(rK^ (7) 

r=0 

If the sequence A is real, then also 

PSD^(5) = VAr^(j)cos^. (8) 

Z / y 

j = 
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4 Complementary sequences 

The associated sequences of the base blocks of an SDS have an important property known as 
complementarity. Let us begin with the general definition of this property. 

Definition 2 Let Ai, . . . ,At be complex sequences of length v. We say that these sequences are 
complementary if 

t 

^PAFa, = [ao, a, ..^.,a ] 
v-i terms 

for some ao and a (the VKF -constants). 

In the next theorem we show that a finite collection of complex sequences has constant sum 
of PAF values at nonzero shifts if and only if the same is true for the PSD values, and we find 
simple formulae expressing the PSD-constants in terms of the PAF- const ants. (This theorem 
is a generalization of ^ Theorem 2].) 

Theorem 2 Let Ai, . . . ,At be complex sequences of length v. These sequences are complemen- 
tary, i.e., 

t 

J]PAFa^ = [ao, (9) 
v-i terms 

if and only if 

t 

^PSD^^ = [/3o, (10) 
v-i terms 

The PAF -constants ao and a and the PSD-constants Po and /3 are related as follows: 

= ao + {v — l)a, /3 = ao — a. (11) 



Proof Note that can be written as: 

t 

PAF(A,) = (ao - a) ■ [1, 0, . . . , 0] + a ■ [1, . . . , 1]. 

i=l 

By applying the DFT and by using the Wiener-Khinchin theorem and the fact that DFT is a 
linear operator, this implies that 

t 

^PSD(Ai) = (ao-a)-DFT[l,0,...,0] + a-DFT[l,...,l] 

j=i 

= (ao - a) ■ [1, . . . , 1] + a ■ [t;,0, . . . ,0] 
= [(3o,l3,...,f3] 
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Hence ( ITOl) holds as well as ( ITTl) . Conversely, assume that ( fTOj) holds. By applying the inverse 
DFT we obtain 

t 

5^PAF(A,) = (/3o-/3)-DFT-i[l,0,...,0]+/3-DFT-i[l,...,l] 

i=l 

^°~^-[l,...,l]+/3[l,0,...,0] 



V 



By solving the equations (ITT]) we obtain that 



□ 



/3o-/3 , ^ /3o-/3 
ao = hp, a = . 

V V 

Note that the PSD values are always nonnegative. Hence, if Ai, . . . , are complementary 
complex sequences of length v with the PSD-constants /3o and /3, then for r = 1, . . . , t we have 

PSD^,(s)</3, s = l,2,...,i;-l. (12) 

We shall refer to this inequality as the PSD-test. 

In the following proposition we show that the associated sequences of the base blocks of an 
SDS are complementary sequences and we compute the PAF and PSD-constants. In a special 
case these constants were computed in [71 Example 3]. 

Proposition 1 Let Ai, . . . ,At be the {±1}- sequences associated to the base blocks Xi, . . . ,Xt 
of an SDS with parameters {v; ki, . . . , kf, A). Then 

t 

N{Ai) = 4n + {tv - 4n)T. (13) 

i=l 

Moreover, the sequences Ai, . . . ,At are complementary with PAF-constants 

ao = tv, a = tv — 4n, (14) 

and PSD-constants 

= tv, 13 = An. (15) 

Proof Note that x^T = T for each i, and so we have = vT and TXi = TX* = kiT. By Eq. 
© we have N{Ai) = (T - 2Xi){T - 2X*) which gives 

N{A,) = {v-4h)T + 4N{X,). (16) 



Summing over alH = 1, . . . ,t and by using Eq. (j4j), we obtain Eq. (jT3|) . Using Eq. ([6]), we 
have 

v-l 

n{a,) = J2p^FaAj)x'. 

j=0 



By summing these equations over i = 1, . . . ,t and by using ( fT3ll . we deduce that 

t 



i=l 



tv, J = 0. 



Hence, the equations ffT^ hold, and by using Theorem [2] we deduce that the equations ffTB]) 
also hold. 

□ 

When searching for SDS, we can discard the trial base block Xr if the associated sequence 
Ar fails the PSD-test: PSD^^(s) > (3 for some s 7^ 0. In this case we can restate the PSD-test 
as follows. 

Lemma 2 Under the hypotheses of Proposition^ for r = 1, . . . ,t we have 



v-1 



y^Nx^{j)cos'^ <n-K, se{l,...,v-l}. (17) 



■ 1 ^ 



Proof Recall that the DFT of a constant function vanishes for all shifts s except for s = 0. 
Hence, by using Eqs. and ([ISD, we obtain that DFT(PAFaJ(s) = 4 ■ DFT(A^x,.)(s). Thus, 
we must have 



v-1 



4 + ^ iVx„(j) COS ^) </3 = 4n, s ^ 0, 



i.e., ([nD holds. 



5 Compression of sequences 

If we have a collection of complementary sequences of length v = dm, then we can compress 
them to obtain complementary sequences of length d. We refer to the ratio v/d = m as the 
compression factor. Here is the precise definition. 

Definition 3 Let A = [ao, ai, . . . , a^_i] be a complex sequence of length v = dm and set 

af^ = aj + aj+d + • • • + aj+(rn-i)d, j = 0,...,d-l. (18) 

Then we say that the sequence A'^^'^ = [a'^\ a^f\ . . . , a^f}^-^] is the m-compression of A. 

Let us now show that the complementarity is preserved by the compression process. At the 
same time we shall compute the PAF and PSD-constants of the compressed sequences. 
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Theorem 3 Let Ai = [aio, an, ... , ai,t,-i], i = 1, . . . ,t, be complementary complex sequences of 
length v = dm with the PAF - constants ao,a- Then the corresponding m-compressed sequences 
A — ['^^o^ • • • ; '^i^d-i]' = 1; • • • ; ^; ^^i"^ Complementary with PAF - constants : 

Oq'^^ = ao + — C(^'^^ = fnci- (19) 

The PSD-constants of the original and compressed sequences are the same. 

Proof By using f lTSj) we compute PAF^{d) (k) for arbitrary shift k: 

i 



PAF^(.)(A;) = 



j=0 

m— 1 ci— 1 



^ ^ ^ ^ ^i,j+k+rd^i,j+sd- 
r,s=0 j=0 



By applying the substitution r — )■ r + s, we obtain 

m— 1 d— 1 



r,s=0 J=0 
m— 1 

= 5^PAFA,(A; + rci) 

r=0 

By summing over z from 1 to t, we obtain 

t m—l t 

^PAF^(.)(fc) = 5^5^PAF^,(fc + rc?). 



r=0 i=l 



Since PAFyi.(A; + rd) is equal to ao if = r = and is equal to a otherwise, we conclude that 
the sequences A^^\ . . . , A^"^^ are complementary and that their PAF-constants are given by Eqs. 
(fT9l) . The assertion about the PSD-constants follows from Theorem O 



□ 

By applying Theorem [3] directly to SDS with parameters {v; ki, . . . , kt] A), we obtain the 
following result. 

Theorem 4 Let Ai, . . . ,At be the sequences associated to the base blocks of an SDS with pa- 
rameters [v] ki, . . . , kf, X) . Suppose that v = dm and let A^f\...,Al-^^ be the corresponding 
m-compressed sequences. Then the PAF-constants of the compressed sequences are given by 

aj,'^^ = m{tv - An) + 4n, a^'^'^ = m{tv - 4n), (20) 

where n is defined by Eq. 
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Proof By Proposition [H the sequences Ai,...,At are complementary with PAF-constants 
given by Eqs. ( TI^ . It remains to apply Theorem [3l 



□ 

Remark 1 Theorem ^ generalizes Theorem 1], which deals with the special case: SDS 
(v; r, s; A) with n = r + s — \ = {v — l)/2. (These SDS are used to construct circulant D-optimal 
matrices.) 

Remark 2 Theorem^ has two useful corollaries for 

• m = 2, in which case {a^fl G {0,2}, 

• m = 3, in which case {a^fl G {1,3}. 

An easy counting argument allows us to compute the number of a^f that assume a specific 
absolute value. 

Corollary 1 With the hypotheses and notation of Theorem^ let m = 2 and denote by vq,V2 
the number of a^j^ terms in all 2-compressed sequences Af''^ that have absolute value equal to 

0, 2 respectively. Then uq = n and 1/2 = td — n. 

Proof Since m = 2, from Eq. (!20l) we have a^^^ = A{td — n). Therefore 0-z/o + 4-i/2 = 4:{td — n), 

1. e., 1/2 = td — n. As uq + 1^2 = td, we have uq = n. 

□ 

Corollary 2 With the hypotheses and notation of Theorem^ let m = 3 and denote by z/i,z/3 
the number of a^f^ terms in the 3-compressed sequences A^f' that have absolute value equal to 
1, 3 respectively. Then Vi = n and = td — n. 

Proof Since m = 3, from Eq. f l2U]) we have a^^^ = 9td — 8n. Therefore 1 ■ z/i + 9 ■ z/3 = 9td — 8n. 

As Ui + = td, we have ui = n and = td — n. 

□ 

6 Computational results for SDS with two base blocks 
6.1 Necklaces, bracelets and charmed bracelets 

When compressing the pair A, B of periodic complementary {±1} sequences of length v = dm, 
with compression factor m, we obtain the pair A^'^\ B^"^^ of periodic complementary sequences 
of length d whose elements are integers belonging to the set {m, m — 2, ... ,2 — m, —m}. The 
first equation of (fT5l) shows that the sum of squares of the elements of A^^^ and B^'''^ together 
is equal to 2v. In our applications we only used the compression factors m = 2 and m = 3. In 
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these two cases the Corollaries [U and [2] determine the total number of squares 0^, 2^ and 1^, 3^, 
respectively. There are several ways to distribute these squares over the sequences A^'^^ and 
B^'^\ and we treat each of them separately. 

We shall use combinatorial objects known as necklaces, bracelets and charmed bracelets 
[TSf [T6] Roughly speaking the necklaces are cyclic arrangements of v objects where we do not 
distinguish the arrangements obtained by cyclic shifts. If we enlarge the equivalence classes by 
allowing also that the arrangement be reversed, then we refer to these equivalence classes as 
bracelets. We can further enlarge these equivalence classes by moving, for each i, the object 
in position i to the position si (mod v) where s is a fixed integer relatively prime to v. We 
shall refer to this operation as a multiplication by s. Note that the multiplication by s = 
— 1 is the same as the reversal. We refer to these enlarged equivalence classes as charmed 
racelets. We emphasize that the multiplication operation does change the PAF of a sequence, 
but when performed simultaneously on a pair of complementary sequences it preserves the 
complementarity property. 

The sequences A and B, as well as the compressed sequences A^'^^ and B^'^\ can be viewed as 
being cyclic arrangements of length v and d, respectively. The cyclic shifts and the reversal do 
not change the periodic autocorrelation of the sequence, and so they can be applied separately 
on A and B without destroying the complementarity property. However, the multiplication by 
s 7^ — 1 (mod v) has to be performed simultaneously, with the same s, on both A and B as 
otherwise the complementarity property may be destroyed. 

In our searches using the compression method we first construct the compressed sequences 
A^'^^ and B^^\ They must have the specified row sums and pass the PSD test. We choose the 
first sequence to be a representative of a charmed bracelet while the second sequence can be 
chosen only as a representative of an ordinary bracelet. Then we proceed to select the pairs 
which are complementary. If there are no such pairs this means that the SDS that we 
are looking for do not exist. Otherwise we examine each of the complementary pairs 
and try to lift them to obtain a complementary pair A, B. 

6.2 The range v < 50 

We shall prove the following theorem resolving the existence of all but one of the 15 open cases 
from the list in the Introduction. The only remaining undecided case is now (49; 21,4; 9). 

Theorem 5 

• There do not exist SDS with following parameters: (41; 15, 6; 6) (43; 9, 4; 2) (44; 19, 2; 8) 
(45; 18, 2; 7) (46; 21, 6; 10) (47; 9, 5; 2) (47; 12, 3; 3) (47; 14, 2; 4) (47; 15, 5; 5) (48; 14, 3; 4) 
(49; 10, 3; 2) (50; 8, 7; 2) (50; 20, 4; 8) 

• There exist SDS with the following parameters: (50; 22, 21; 18) 
Proof We proceed by examining the 14 {v; r, s; A) cases one by one. 
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1. For the case (41; 15, 6; 6) the apphcation of the PSD test with the constant value 60 gave 
1040 normahzed candidate A-sequences and 13104 normahzed candidate B-sequences. 
Subsequently, there was no match found between these two sets of candidate sequences. 

2. For the case (43; 9, 4; 2) there are no A-sequences that pass the PSD test with the constant 
value 44, therefore we conclude immediately that such SDS do not exist. 

3. For the case (44; 19, 2; 8) there are no A-sequences that pass the PSD test with the 
constant value 52, therefore we conclude immediately that such SDS do not exist. 

4. (45; 18, 2; 7) there are no A-sequences that pass the PSD test with the constant value 52, 
therefore we conclude immediately that such SDS do not exist. 

5. For the case (46; 21, 6; 10) we used 2-compression to prove that there do not exist such 
SDS. Here are the details of this computation. First note that in this case we have 
n = r + s — A = 21 + 6 — 10 = 17. Taking m = 2, c/ = 23 in Theorem HI we see that we 
need to find two sequences of length 23 with elements from {—2,0, +2} such that their 
PSD values add up to 4n = 4 ■ 17 = 68. Since we know from corollary [1] that the total 
number of elements in the two sequences is equal to uq = 17 and that the total number 
of ±2 elements in the two sequences is equal to z/2 = 29, we deduce that there are four 
cases to consider, using ordinary and charmed bracelets [U]. For each one of these four 
cases: (1) we compute the charmed bracelets that give 2-compressed sequences of length 
23 all of whose PSD values are smaller than 68 and (2) we compute the ordinary bracelets 
that give 2-compressed sequences of length 23 all of whose PSD values are smaller than 
68. We summarize the results in the following tables. 



A-sequences 



Case 


charmed bracelets 


charmed bracelets passing PSD 


1 


2,116,296 


85 


2 


475,020 


2,009 


3 


54,264 


4,552 


4 


3,015 


1,442 



B-sequences 



Case 


ordinary bracelets 


ordinary bracelets passing PSD 


1 


2,277 


1,749 


2 


3,685 


1,419 


3 


1,210 


22 


4 


44 






Notice that case 4 can be eliminated at this stage, since there are no ordinary bracelets 
that pass the PSD test. The next step is for each of the remaining three cases to look for 
pairs of sequences that have constant PAF and we summarize the results in the following 
table (where the symbol — > indicates removal of sequences with duplicate PAF) 
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Case 


A- sequences 


B-sequences 


# of pairs 


1 


85 ^ 84 


1,749 ^ 1,716 


39 


2 


2,009 ^ 1,970 


1,419 1,419 


34 


3 


4,552 4:,m 


22 22 






Notice that case 3 can be eliminated at this stage, since there are no pairs of matching 
sequences at all. Subsequently wc used the 39 + 34 = 73 matching pairs of sequences 
of length 23 with elements from {—2,0, +2} to construct all the corresponding pairs of 
uncompressed sequences of length 46 and check whether they form an SDS with the 
required parameters. We did not find any such SDS. 

6. For the case (47; 9, 5; 2) there are no A-sequences that pass the PSD test with the constant 
value 48, therefore we conclude immediately that such SDS do not exist. 

7. For the case (47; 12, 3; 3) there are no A-sequences that pass the PSD test with the 
constant value 48, therefore we conclude immediately that such SDS do not exist. 

8. For the case (47; 14, 2; 4) there are no A-sequences that pass the PSD test with the 
constant value 48, therefore we conclude immediately that such SDS do not exist. 

9. For the case (47; 15, 5; 5) there are no A-sequences that pass the PSD test with the 
constant value 60, therefore we conclude immediately that such SDS do not exist. 

10. For the case (48; 14, 3; 4) there are no A-sequences that pass the PSD test with the 
constant value 52, therefore we conclude immediately that such SDS do not exist. 

11. For the case (49; 10, 3; 2) there are no A-sequences that pass the PSD test with the 
constant value 44, therefore we conclude immediately that such SDS do not exist. 

12. For the case (50; 8, 7; 2) the application of the PSD test with the constant value 52 gave 
1130 normalized candidate A-sequences and 2910 normalized candidate B-sequences. Sub- 
sequently, there was no match found between these two sets of candidate sequences. 

13. For the case (50; 20, 4; 8) we used 2-compression to prove that there do not exist such 
SDS. 

14. We give four non-equivalent examples of SDS (50; 22, 21; 18): 

{0, 1, 2, 3, 6, 7, 9, 13, 14, 16, 18, 20, 22, 23, 26, 27, 30, 35, 37, 41, 45, 46} 
{0, 1, 2, 3, 4, 5, 6, 8, 11, 12, 14, 17, 20, 22, 29, 30, 32, 37, 38, 39, 42} 

{0, 1, 2, 3, 4, 6, 7, 8, 9, 14, 16, 18, 20, 21, 25, 31, 32, 35, 36, 42, 44, 45} 
{0, 1, 2, 4, 5, 8, 9, 10, 12, 14, 18, 21, 23, 24, 27, 29, 32, 34, 35, 39, 42} 

{0, 1, 2, 3, 5, 8, 9, 11, 14, 15, 19, 21, 24, 25, 29, 30, 32, 36, 38, 39, 41, 43} 
{0, 1, 3, 5, 6, 7, 8, 9, 10, 13, 16, 18, 20, 21, 24, 25, 31, 32, 33, 37, 41} 
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{0, 2, 3, 4, 6, 9, 10, 12, 13, 17, 19, 20, 24, 25, 28, 29, 30, 33, 38, 39, 41, 47} 
{0, 1, 3, 5, 6, 7, 8, 10, 12, 13, 14, 17, 20, 22, 24, 28, 32, 37, 38, 39, 40} 

All four examples are in the canonical form defined in |5] and since they are different, 
this implies that they are non-equivalent. Note that since 2u — 4n = 2 ■ 50 — 4 • 25 = 0, 
these SDS give rise to {±1} sequences of length 50 with PAF zero, because of ( |T4l) . 

□ 

6.3 The range v > 50 

In this section we consider the SDS (f ; r, s; A) with v = 2n, where n = r + s — X. In particular v 
is even. Recall from Proposition [T] that these SDS give rise to binary periodic complementary 
pairs, i.e., pairs of {±1} sequences of length v with PAF constant a = 0, because of ffH]) . 

Since 52 is a Golay number, there exist Golay pairs of length 52. At the same time they 
are also binary periodic complementary pairs. Since the Diophantine equations + = 2 ■ 54 
and + y'^ = 2 ■ 56 do not have solutions, there are no binary periodic complementary pairs of 
lengths 54 and 56. Therefore, the first interesting value of f > 50 is f = 58, since this value is 
not excluded by the restriction of Arasu and Xiang [T]. In this paper we find 4 non-equivalent 
examples of SDS (58; 27, 24; 22), for the first time. These SDS give 4 respective pairs of binary 
sequences of length 58 with zero periodic autocorrelation function. 

We give four non-equivalent examples of SDS (58; 27, 24; 22): 

{0, 1, 2, 3, 4, 7, 8, 10, 11, 12, 13, 16, 18, 20, 24, 26, 29, 31, 32, 33, 36, 38, 43, 46, 47, 50, 53} 
{0, 1, 2, 3, 7, 8, 10, 11, 12, 13, 16, 17, 21, 22, 24, 27, 30, 34, 41, 42, 43, 45, 47, 49} 

{0, 1, 2, 3, 5, 6, 7, 9, 11, 12, 14, 15, 17, 19, 23, 24, 25, 26, 29, 32, 33, 39, 40, 43, 45, 48, 52} 
{0, 1, 2, 3, 4, 5, 9, 11, 14, 15, 16, 18, 22, 26, 27, 31, 32, 34, 37, 39, 41, 42, 45, 51} 

{0, 1, 2, 3, 5, 8, 9, 11, 12, 13, 14, 18, 19, 21, 24, 25, 27, 29, 32, 34, 35, 39, 41, 43, 44, 48, 49} 
{0, 2, 3, 4, 6, 8, 10, 13, 16, 17, 19, 20, 21, 25, 28, 29, 32, 33, 34, 39, 40, 41, 43, 46} 

{0, 2, 3, 4, 6, 7, 8, 10, 11, 14, 16, 17, 18, 20, 23, 25, 26, 28, 31, 32, 36, 37, 38, 41, 42, 47, 49} 
{0, 1, 2, 3, 5, 8, 9, 10, 12, 16, 17, 18, 22, 25, 28, 30, 35, 37, 41, 44, 45, 46, 48, 49} 

All four examples are in the canonical form defined in [5] and since they are different, this 
implies that they are non-equivalent. 

7 Acknowledgements 

The authors thank Joe Sawada and Daniel Recoskie for sharing improved versions of the their C 
code for computing ordinary and charmed bracelets. The authors wish to acknowledge generous 
support by NSERC. This work was made possible by the facilities of the Shared Hierarchical 
Academic Research Computing Network (SHARCNET) and Compute/Calcul Canada. 



13 



References 



[1] K. T. Arasu and Q. Xiang, On the existence of periodic complementary binary sequences 
Des. Codes Cryptogr. 2 (1992), 257-262. 

[2] L. D. Baumert, Cyclic difference sets. Lecture Notes in Mathematics, Vol. 182 Springer- 
Verlag, Berhn-New York 1971. 

[3] J. H. E. Cohn, A D-optimal design of order 102. Discrete Math. 102 (1992), 61-65, 

[4] D. Z. Dokovic, Note on periodic complementary sets of binary sequences, Des. Codes 
Cryptogr. 13 (1998), 251-256. 

[5] D. Z. Dokovic, Cychc {v; r, s; A) difference famihes with two base blocks and v < 50. Ann. 
Comb. 15 (2011), 233-254. 

[6] D. Z. Dokovic, I. S. Kotsireas, New results on D-optimal Matrices. J. Combin. Designs, 20 
(2012), 278-289. 

[7] R. J. Fletcher, M. Gysin and J. Seberry, Application of the discrete Fourier transform to 
the search for generalised Legendre pairs and Hadamard matrices. Australas. J. Combin. 
23 (2001), 75-86. 

[8] D. Jungnickel, A. Pott, K. W. Smith, Difference sets, in Handbook of combinatorial de- 
signs, 2nd ed. C. J. Colbourn and J. H. Dinitz (eds) pp. 419-435. Discrete Mathematics 
and its Apphcations (Boca Raton). Chapman & Hall/CRC, Boca Raton, FL, 2007. 

[9] H. Kharaghani and W. Orrick, D-optimal matrices, in Handbook of Combinatorial Designs, 
2nd ed. C. J. Colbourn, J. H. Dinitz (eds) pp. 296-298. Discrete Mathematics and its 
Apphcations (Boca Raton). Chapman & Hall/CRC, Boca Raton, FL, 2007. 

[10] I. S. Kotsireas, and C. Koukouvinos, Periodic complementary binary sequences of length 
50. Int. J. Appl. Math. 21 (2008), 509-514. 

[11] S. Kounias, C. Koukouvinos, N. Nikolaou and A. Kakos, The nonequivalent circulant D- 
optimal designs for n = 2 mod 4,n < 54,n = 66. J. Combin. Theory Ser. A 65 (1994), 
26-38. 

[12] L. Martinez, D. Z. Dokovic, A. Vera-Lopez, Existence question for difference families and 
construction of some new families. J. Combin. Designs, 12 (2004), 256-270. 

[13] R. Mathon and A. Rosa, 2-{v, k, A) Designs of Small order, in Handbook of Combinatorial 
Designs, 2nd ed. C. J. Colbourn, J. H. Dinitz (eds) pp.. 25-58. Discrete Mathematics and 
its Applications (Boca Raton). Chapman & Hall/CRC, Boca Raton, FL, 2007. 

[14] D. W. Ricker, Echo Signal Processing. Series: The Springer International Series in Engi- 
neering and Computer Science, Vol. 725 Springer, 2003. 



14 



[15] J. Sawada, Generating bracelets in constant amortized time. SIAM J. Comput. 31 (2001), 
259-268. 

[16] J. Sawada, A fast algorithm to generate necklaces with fixed content. Theoret. Comput. 
Sci. 301 (2003), no. 1-3, 477-489. 

[17] D. R. Stinson, Combinatorial designs. Constructions and analysis. Springer- Verlag, New 
York, 2004. 

[18] Yang, C. H. On Hadamard matrices constructible by circulant submatrices. Math. Comp. 
25 (1971), 181-186. 



15 



